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1.  Introduction 


A  sample  ...  ,  x^"^  is  drawn  from  the  normal  distribution 

(1)  (27  (2)  (2) 

N(p  ,  E),  and  a  sample  x  ,  . .  .  ,  x__  is  drawn  from  N(p  ,  E). 

(1)  (2) 

The  p-component  mean  vectors  p  and  p  and  the  common  covariance 

matrix  E  are  unknown;  it  is  assumed  that  p^^  ^  p  '2^  and  E  is  non¬ 
singular.  Another  observation  x  is  drawn.  It  is  desired  to  classify 
this  observation  as  coming  from  N(p^,  E)  or  N(p^2\  E).  [See  T.  W. 
Anderson  (1951)  or  T.  W.  Anderson  (1958),  Chapter  6.J 

The  observation  x  may  be  classified  by  means  of  the  classification 
statistic 


(1) 

(2) 


W-  -;«>>•  s-1  lx  -  i  ff(1>  +  *«>):  , 

~  <v  <v  ^  4  ~  ~ 


N  N 

-(1)  1  y  (1)  -(2)  1  y  (2) 

x  =  r~  )  x.  ,  X  =  —  )  x. 

N.,  ,u-  -n  ~  N„  ,L.  ~i 

1  3=1  J  2  3=1  J 


(3)  nS  -  I'  (x^-x^)  (xf1)-x(1))'  +  }  (x<2>-x<2>)  <x< 2>-S?<2>)'  , 

~  j=1  ~3  ~  -J  -  j=1  ~  ~ 

and  n  =  N^+N^  -  2.  The  rule  is  to  classify  x  as  coming  from 

N  (p  ^ ,  E)  if  W  >  c  and  from  N(p^\  E)  if  N  _<  c  ,  where  c  may 

be  a  constant,  particularly  0,  or  a  function  of  x^\  x^2\  and  S. 

(1)  (2) 

The  distribution  of  W  depends  on  the  parameters  p  ,  p  , 
and  E  through  the  squared  Mahalanobis  distance 


(4) 


a  =  (p(1)  -  p(2))«  E"1  (p(1)  -  p(2))  , 


which  can  be  estimated  by 


(5) 


a.  (x(1)  -  ;(2V  S'1  <x(1)  -  x(2)) 
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The  limiting  distribution  of  W  as  00  and  ^  00  is  normal  with 

variance  a  and  mean  ^  a  if  x  is  from  N(]i^^,  E)  and  mean  -  \  a 

Z  ~  ~  ~  z 

(2) 

if  x  is  from  N(]iv  ,  E).  Bowker  and  Sitgreaves*  (1961)  for 
and  Okamoto  (1963)  [with  correction,  Okamoto  (1968)]  gave  asymptotic 
expansions  of  the  distribution  of  (W  -  a)/Va  for  x  coming  from 
N(]i^\  E)  and  (W  +  -j  a)/v^a  for  x  coming  from  N(y^\  E)  to  terms 
or  order  1/N^,  l/N^,  and  l/n^  when  ->  00 ,  ■>  °°,  and 

a  finite  positive  constant.  In  particular,  Pr{w  _<  0}  was  evaluated. 

The  statistician,  who  wants  to  classify  x,  may  take  c  to  be  a 
constant,  perhaps  0,  and  accept  the  pair  of  misclassif ication  probabilities 
that  result.  The  asymptotic  expansion  of  the  distribution  of  (W  +  a)/Va 
gives  approximate  evaluations  of  these  probabilities,  which  are  functions 
of  the  unknown  parameter  a  as  well  as  of  c. 

On  the  other  hand  the  statistician  may  want  to  determine  the  cut-off 
point  c  to  adjust  the  probabilities  of  misclassif ication.  Since  the 
limiting  distribution  of  ( and  (W+oO/v'oT  are  N(0,  1)  when 
^x  =  and  £  x  =  ]i  ^ ,  respectively,  a  first  approximation  to  the 

pair  of  misclassif  ication  probabilities  is  $(y  a  +  cv'a)  and  $(—  7^  a  +  cv^t), 
where  $(a)  is  the  cumulative  distribution  function  of  the  standard 
normal  variate.  Since  a  is  an  estimate  of  a,  one  might  base  his 
choice  of  c  on  the  fact  that  the  limiting  distribution  of  (W  -  y  a)//a 
and  (W  +  i  a)//a  are  N(0,  1)  when  ^x  =  and  ^,x  = 

respectively.  In  this  paper  we  make  asymptotic  expansions  of  the  distri¬ 
bution  of  (W  -  y  a)/v^a  and  (W  +  y  a)//a  in  these  two  cases,  respectively. 

and  a^2  should  be  replaced  by  -a^  an<^  -a32’ 


*The  coefficients  a^ 
respectively. 
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2.  The  Asymptotic  Expansion 

The  statistics  x,  x^,  x^2\  and  S  are  independently  distri¬ 
buted  according  to  N(y,  E),  N[y^2\  (1/N^)E],  N[y^2\  (l/N^)^],  and 

W(E,  n) ,  respectively;  here  y  =  %,x  and  W(E,  n)  denotes  the  Wishart 
distribution  with  n  degrees  of  freedom.  We  write 


(6) 

Then 


W 


-  -  a  =  (x(1)  -x<2>V  s"1 


x^O'  S  (x  -  xVJ-y)  . 


(1). 


(7)  pr. 


W-|a 


f 


<  u)  =  Pr  ;  (x(1)-x^2))'  s"1  (x-lij 


<  u  /<xa>-x<2V s-rS(1>-i(2))  + 

—  V  ~  ~  ~  ~  ~  ~  ~  ~  ~  ~  / 


Since  x  has  the  distribution  N(y,  E)  independently  of  x^,  x^ 
and  S,  the  conditional  distribution  of  (x ^-x ^ ) ' S  ^(x-y)  is 


N[0,  (x(1)-x(2)),S_1Zs  1(x(1)-x(2))],  and 


(8) 


(x^-x^)  's  1(x-y) 


r  = 


/fe(1)-x(2))'S-1ES“1S(1)-i(2)) 


has  the  distribution  N(0,  1).  Then  (7)  is 


W  -  j  a 

(9)  Pr  ( -  <  u^>  =  Pr  (  r  < 


L 


& 


u/Cx^-x^)  *S  1(x^1^-x^2))+(x^1^-x^2^)'s_1(x<'1^-y) 
/(x^^-x^^S^ZS^Cx^^-x^2^) 


=  £$ 


u/u<1)^(2))*s"1fr<1U<2))  +  (x(1)-x(2))'S-1(x(1)-y) 


/(x^-x^)  ,S-1ZS_1(x^1^-x^2^ 


where  the  expectation  is  with  respect  to  x^\  x^2\  and  S. 
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The  distribution  of  W  and  a  is  invariant  with  respect  to  the 

and 


*  *(1)  (1) 

transformations  x  =1  A  x  +  b,  x.  =  A  x.  +  b,  j  =  1,  . . .  ,  N1  , 


-J  ~ 

* (2)  (2) 

x.  =  A  x.  +  b  ,  where  A  is  nonsingular.  The  maximal 

J  ~  ~ 

invariant  of  these  transformations  is  the  distance  a,  given  by  (4).  We 

can  choose  A  and  b  to  transform  E  to  I,  to  6  =  (A,  0 . 0), 

where  A  =  v'ct,  and  y,  to  0.  We  shall  first  treat  the  case  where  y  =  y.  . 

~  X  ~  ~  ~  X 

The  vectors  x^^  and  x^^  and  the  matrix  S  are  distributed 

—Q's  _ (2)  —(1)  ' 

independently.  The  distribution  of  [  (x  -x  )\  x  '  ]’  is 


(10) 


N 


*6  ' 


k + k> 

—  i 
Ni~ 


Let  Y,  Z  and  V  be  defined  by 


k -1' 


N, 


I 

J  • 


(11) 


x(D_x(2)  =  6 ,  +  —  y  ,  x(1) 

/n 


/z 


(12) 


s  =  I  +  —  v  . 
■  y£  ~ 


Then  the  joint  distribution  of  (Y1 ,  Z’)1  is 


(13) 


N 


\9/ 


(1  .  1  N 
nk + ^ 


\ 


n 

N., 


n  A 


N  5 


n 

N, 
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Then  (9)  is 


W  -  f  a 

(14)  Pr  < - - —  <_  u 

Ja 


= 


uA6  +  —  Y) '  (I  +  —  V)  1  (6  +  —  Y)  +  —(6  +  —  Y)  '  (I  +  —  V)-1Z 


/n 


& 


/n  ~  )/ri  ~  /n  -/n 


/(6  +  —  Y)  ’  (I  +  —  V)“2(6  +  —  Y) 


& 


& 


We  can  write 


(15)  (I  +  —  V)  1 

~  /S  ~ 


I  —  ~  V  +  “  V2  —  —pr  v3l^v4-  --JJ  v5  (I  +  —  V)-1  , 

-  r-  -  n  ~  3/2  -  2  ~  5/2  ^  — 

/n  n  n  n 


(16)  (I  +  —  V)  2  =  I~—  V+-V2-  v3  +  -~r  v4 
~  ^  ~  ~  &  ~  n  ~  n3/2  ~  n3/2  ~ 


n 


k(6V5+  -  V6)  (I  +  —  V)"1  . 
5/2  ~  ~  ~  y£  ~ 


Then  (as  Taylor  series  expansions)  we  have 


(17) 


(6  +  —  Y)’  (I  +  —  V)_1(<5  +  —  Y) 


/n  /n 


/n 


1/2 


6 ' 6  +—  (26’Y-6’V6)  +—  (S’V2S  +  Y'Y  -  26'VY)  +  r.  (Y,Z,V) 

~  ~  r~  ~  ~  ~  — -  n  IV  Arf  A/  Arf  I  T~l  #V  A/  IV 

/n 


1/2 


=  A  +  — i —  (26 'Y  -  6'V6)  +  — 
2A  /n . n 


(6’V26  +  Y’Y  -  26’VY) 


(26'Y  -  6’V6)2 
8A  ~  ~  ~  ~~ 


+  r2n(I’?’Y>  ’ 
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<  < 


(18)  —  (5  +  —  Y)’  (I  +  —  V)-1  Z  =  —  6’Z  +  -  (Y’Z  -  6'VZ)  +  r0  (Y,Z,V)  , 

7  ~  T  '  ~  'v  /  *v»  *v>  T  *V  <v  n  >V  <V  'W'v  'V/  >v  <v 

vn  /n  yn  /n 


(19)  (6  H - —  Y) ’ (I  +  —  V)  2  (6  +  —  Y)1 

•V  7"  *  -V  <v  /  »v  <v  /  >v  I 

^  i/r»  t/ r> 


=  [S'S  +—  (26 ’Y-26 ’V6)  +  -  (36’V26  +  Y’Y  -  46’VY)  +  r.  (Y,Z,V)1 

'V  /  *v  'V  »V  *  '  n  'V  »V  #V  «V/  *v/  Ai  ifn  AJ  Ay  A/  J 


A  1/2 


=  -  -  — —  (6  ’Y-6  ’V6)  -  -  f-^-  (36’v26  +  Y’Y  -  46’VY) 

A  «  -J  /  —  A/  a/  Ay  Ay  Ay  XT  !  a  A  J  Ay  Ay  Ay  Ay  Ay  Ay  Ay  Ay 

A  /n  ]_2A 


-  (6’Y-6'V6)2  +  r,  (Y,Z,V) 

1  Ay  Ay  *vy  *vy  Ay  J  Tl  Ay  Ay  Ay 

2A  J 


Here  rjn(Y,Z,V),  j  =  1,  ...  ,  5  is  a  remainder  term  consisting  of  1/n  ' 

times  a  homogeneous  polynomial  (not  depending  on  n)  of  degree  3  in  the 

2 

elements  of  Y,  Z,  and  V  plus  1/n  times  a  homogeneous  polynomial  of 

-5 /2 

degree  4  plus  a  remainder  term  which  is  0(n  )  for  fixed  Y,  Z,  and 


The  argument  of  $(  )  in  (14)  is  the  product  of 


(20)  uA  +  —  -rr[  (26  ’Y-6  ’V6)+6’Z]  +  -  [77  (6  ’V26+Y’Y-26  ’VY) 

/  yC /\  Ay  Ay  Ay  a-  Ay  Ay  XT  Ay  Ay  Ay  Ay  Ay  Ay  A/Ay 

vn 


(26’Y-6’V5)2  +  Y’Z-S’VZ]  +  r,  (Y,Z,V) 

_  ,  J  Ay  *vy  Ay  Ay  Ay  Ay  Ay  Ay  AyAy  0X1  Ay  Ay  Ay 

oZ\ 


and  (19),  which  is 


u  +  —  (-—  6’V6  +  \  6’Z)  +  -  [^r  (6’Vy-6’V26)  +  ^7  (~6’Y6’V6  +  \  (6’V6)2) 


v'n  2A 


n  V 


+  —  Y’Z  -  \  6’VZ  -  ir-  6’YZ’6  +  \  6’Z6’V6]  +  r7  (Y,Z,V) 

/  \  Ay  Ay  jj[  Ay  Ay  Ay  ^  J  Ay  Ay  Ay  Ay  ^  J  Ay  Ay  Ay  Ay  Ay  /XT  Ay  Ay  Ay 


=  u  +  —  C(Z,V)  +  -  D(Y,Z,V)  +  r7  (Y,Z,V)  , 

~  ~  IT  Ay  Ay  Ay  /  XT  Ay  Ay  Ay 


6 


say  [as  the  definition  of  C(Z,V)  and  D(Y,Z,V)]  and 
r^n(Y,Z,V)  have  the  same  properties  as  r^n(Y,Z,V),  j 


r,  (Y,Z,V)  and 
on  ~  ~  ~ 

“  1  j  •••  9  3a 


A  Taylor  series  expansion  of  $(  )  in  (14)  gives 


*[u  +  —  C(Z,V)  +  i  D (Y, Z ,V)  +  r7  (Y,Z,V)] 
j—  ~  ~  n  ~  ~  ~  /n~~~ 

vn 


(22) 


=  #(u)  +  <|>(u)  C(Z,V)  +  -  [D (Y,Z,V)  -  in  C2(Z,V)]J 

/£  ~  ~  n  -  2 

+  "3/2  r8(~’?’~)  +i2  r9^,Z,V)  +  r10n(Y,Z9V)  , 
n  n 

where  rg(Y,Z,V)  is  a  homogeneous  polynomial  (not  depending  on  n  but 
depending  on  u)  of  degree  3  in  the  elements  of  Y,  Z,  and  V,  r^(Y,Z,V) 
is  a  polynomial  of  degree  4,  and  r.^n(Y,Z,V)  is  a  remainder  term,  which 
is  0(n  ^2)  for  fixed  Y,  Z,  and  V  (and  u) . 

Let  be  the  set  of  Y,  Z,  and  V  such  that  | y^ |  <  2/log  n, 


\z±\  <  2/log  n,  i  =  1,  ...  ,  p  ,  and  | v±^ ]  <  2  log  n,  i,  j  =  1,  ...  ,  p. 
As  shown  in  the  Appendix, 


(23) 


Pr{J  }  =  1  -  o(n  )  . 


n 


The  difference  between  £$(  )  and  the  integral  of  $(  )  times  the  density 

_2 

of  Y,  Z,  and  V  over  J  is  o(n  ),  because  0  <  $(  )  <  1.  In  J 
~  ~  ~  n  —  —  n 

each  element  of  Y,  Z,  and  V  divided  by  /n  is  less  than  a  constant 


times  (log  n//n)  .  Hence 
(24) 


rl0n(Y,Z,V)  <  constant  x  ^ 

/n 


-2 

and  the  integral  of  this  times  tne  density  of  Y,  Z,  and  V  over  is  o(n  ), 
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Since  fourth-order  absolute  moments  of  Y,  Z,  and  V  exist  and  are 
bounded,  the  integral  of  rQ(Y,Z,V)  times  the  density  of  Y,  Z,  and 

V  over  is  bounded;  hence,  the  contribution  of  this  term  (with  the 

-2  -2 
factor  n  )  is  0(n  ). 

-1 /?  -1  12 
The  differences  between  n  ^C(Z,V),  n  [£  D(Y,Z,V)  -  —  uC  (Z,V)] 

—3/2  —1/2 

and  n  rg(Y,Z,V)  and  the  integrals  over  J  of  n  /  C(Z,V), 

—1  12  -3/2 

n  [D(Y,Z,V)  -  y  uC  (Z,V)]  and  n  rg(Y,Z,V)  times  the  density  of 

_2 

Y,  Z,  and  V,  respectively,  are  0(n  ).  Thus 


=  $(u)  +  cj>(u) 


£C(Z,V)  [£D(Y,Z,V)  -|  £C2(Z,V)] 


) 


+  ^3/2  &  +  0(n_2> 


=  $(u)  +  <j>(u)  {—  $  C(Z,V)  +  ~  I  &D(Y,Z,V) 

J n 


-  |  ^C2(Z,V)3>  +  0(n"2)  . 


because  the  third-order  moments  of  the  elements  of  Y,  Z,  and  V  are 
either  0  or  0(n  2). 

Since  C(Z,V)  is  linear  and  homogeneous ,  j£C(Z,V)  =  0.  Since 
(Y,Z)  and  V  are  independent 


(26)  £D(Y,Z,V)  =  -  ~  l  6'V26  +  j  ~  £(<5’V6)2  +  j  £Y’Z  -  —  £6'YZ'<5 

A  A 


u  t  \  ,  7  u  „.4  ,  1  n 

72  A  (P+1)  +  8  TV  2A  +  A  N~  p 


-(p  -  |)u  +  (P-1)  f-i 


i-  A2  S_ 
3  N 
AJ  W1 
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since 


(27)  £S'V26  =l6,VV,6  =  A2  t  Jj  vli  A  V11  +  .£  ?jVli)  A  (p+1)  ’ 

~  ~  ~  ~  ~~  ~  1  _L  1—  Z 


(28) 


£,(S’VS)2  =  A4  &v2  =  2A4  . 
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We  have 


(29) 


£c2(Z,V)  =  — -  ^(6’V6)2  +  £S*ZZ'6 


4A 


Thus 


u  OA4  1  n  a2  1  2  n 

4*2A  +  2N  A  ~  2  U  +N 

4A  AZ  N1  ^  W1 


(30)  £D(Y,Z,V)  -  §  ?CZ(Z,Y)  -  (p-1)  §-  |  -  (p  -  f  +  f  i  )u  -  i  Uj  . 


3  .  1  n 


1  3 


Replacing  n/N^  by  its  limit  1+k,  we  have 


(31)  Pr  {  =  $(u)  + -~  <KU) 


/a 


(p-1)  v  ,  1  ,  1  .  v  13 

-  (1+k)  -  (P  “  T  +  Y  k)u  -  u 

l  /a 


+  0 (n  2) 


when  A  x  =  y 


(1) 


Interchanging  and  gives 


(32)  Pr 


w  +  ia 

/a 


£  v)  =  <j)(v)  -  —  4> (v) 


P-1  „  .  U  .  ,  1  .  1  v  ,1  3 

r  (1  +  k}  +  (P  -  4  +  2k)v  +  4  V 

,/a 


+  0(n"2)  , 


when  £  x  =  u^2^ 
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3.  Discussion 


If  and  costs  of  misclassif ication  are  equal,  the  minimax 

classification  procedure  is  defined  by  the  cut-off  point  0  for  W;  a 
cut-off  point  different  from  0  will  increase  one  probability  of  mis¬ 
classif  ication  and  decrease  the  other.  The  inequality  W  _<  0  is 

1/—  1?—  It-  i  _ 

equivalent  to  (W  -  a) //a  ^  -  —  va  ,  and  -  —  /a  estimates  -  —  /a 

=  -  y  A.  For  most  purposes,  then,  one  is  interested  in  u  _<  0.  Then 
the  correction  term  to  <&(u)  is  nonpositive;  use  of  the  normal  approxi¬ 
mation  alone  tends  to  underestimate  the  probability  of  misclassif ication. 

The  correction  term  decreases  as  the  distance  A  between  the  two  popula¬ 
tions  increases  if  p  >  1  and  for  nonpositive  u  the  correction  term 
increases  with  the  number  of  coordinates  p  (for  fixed  A) . 

The  expansions  of  Pr{(W  -  y  <_  u|y=y^^}  and  Pr{(W  +  — 

u]y=y'  s  given  by  Okamoto  (1963)  can  be  obtained  by  the  method  of 
this  paper.  It  is  interesting  that  the  expansions  for  (W  + -|-  i^a) //a  here 
are  much  simpler  than  the  expansions  for  (W  +  y  v/ct)/Va  as  given  by 
Okamoto.  At  ]i  =  --yA  =  -yv/,a  (corresponding  to  the  cut-off  point  0) 
the  correction  term  of  order  .1/n  to  the  probability  for  (W  +  y 
is  about  y  as  much  as  for  (W  +  y  fa./ /a . 

As  indicated  in  the  introduction,  the  statistician  may  want  to  use  the 
evaluation  of  Pr{(W  -  y  /a) //a  u}  in  order  to  set  the  cut-off 

point  c  =  uv'a”  +  y  /a  in  order  to  obtain  a  specified  probability  of 
misclassif ication  or  at  least  approximate  a  specified  probability.  The 
crudest  approximation  is  to  take  u  so  $(u)  is  the  specified  probability. 
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This  approximation,  however,  is  not  very  good;  the  error  of  the  approxi- 

3/2 

mation  is  evaluated  above  to  order  1/n  .  The  error  depends  on  the  unknown 

parameter  if  p  >  1.  To  get  a  better  approximation  let  $(u+Au)  be  the 
specified  probability,  where 

a  1  (p-l)(l+k)  ,  1.1,,  13 

Au  - — — - -  -  (p  -  7-  +  —  k)u  -  7-  u 

n  ^  4  2  4 

Then  the  actual  probability  is  the  specified  one  with  an  error  of  order 

~3/2  -1/2 

n  (because  /a  is  with  an  error  of  order  n  ) . 

For  further  discussion,  see  Anderson  (1972). 


Acknowledgements .  The  author  is  indebted  to  John  Miller  and  Fred 
Nold  for  assistance  in  preparing  this  paper. 
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APPENDIX 


To  control  the  errors  of  approximation ve  define  the  set  by 

|Yj|  <  2 /log  n,  I  z  I  <  2/log  n,  j  =  1 ,  p ,  and  |v_|  <  2  log  n, 

i,  i  =  1,  ...  ,  p.  We  want  to  show  Pr{j  }  =  1  -  0(n  ). 

n 

We  have  when  E  =  I 

_  1 

(A.l)  Pr{|y.l  >  2/log  n)  =  — —  f  e  2  dv 

^  /2Tr  2/log  n 

-  j  (2/log  n)2 

<  - - - 

/2tt  log  n 

=  1 

n2  /2tt  log  n 

=  °(n  2) 


by  use  of  Mill's  ratio.  Then 
(A. 2)  Pr{|yj|  <  2/log  n, 


Now  consider  V 
when  2  =  I  is 


(V‘ 


|z_.  |  <2  log  n,  j  =  1,  ...  ,  p} 

=  1  -  o(n  2)  . 

The  moment  generating  function  of 


nS 


(A.  3) 


%> 


tr0nS  nT .  ?  . 

>  "  =  u  1°'1 


9.  .s 

ij  ij 


=  | I  -  20 


-  1/2  n 


where  0=0'.  We  use  the  Tchebychef f-type  inequality  [Chernoff  (1952),  for 
example]  for  an  arbitrary  random  variable  X  and  0  >  0 
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(A. 4) 


■0  a 


>  Pr{X  >  a}  . 


Then 

(A.  5) 

for  0  <  0 

(A. 6) 

Similarly 
(A.  7) 

for  0  <  0 

(A.  8) 

Similarly 
(A. 9) 


e 


/o  e 


0 (X-a ) 


Pr{v.  .  >  2  log  n}  =  Pr{v/n  s .  .  -  /n  >  2  log  n} 
xi  ii 

=  Pr{n  s..  >  n  +  2v^n  log  n) 

ii  ° 

<  (l-20)-^ ^  n  e-®(n+2v/^  log  n) 


Let 


0  =  k/y'ri,  where  k  >  1. 


For 


n  >  4k 


2 


Pr{vJ4  >  2  log  n>  <  (l-2k/Vn)-1//2  n  e  2k  log  n-k^ 

-2k  log  n 


ii 


<  constant  x  e 


=  0(n“2k) 
=  o(n  2). 


-2 

Pr{-v  ^  >  2  log  n}  =  o(n  ).  We  have  for  i  j 


Pr{v..  >  2  log  n}  =  Pr{v^T  s .  ,  >  2  log  n} 
ij  ij 

=  Pr{n  s.  .  >  2v/ri  log  n} 
ij 

<  e-02/n  log  n  (1_02-j“l/2  n 


1  2 
<  — .  Let  0  =  k/Vr i,  where  k  >  1.  For  n  >  4k 


Pr{v^  >  2  log  n}  <_  constant  x  e 

=  o (n  2). 

-2 

Pr{-v„  >  2  log  n)  =  o(n  ).  Then 


-2k  log  n 


Prt|v  |  <  2  log  n,  i,  j  =  1,  ...  ,  p}  =  l-o(n2). 
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